Strongly damped nuclear collisions: zero or first sound ? 
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The relaxation of the collective quadrupole motion in the initial stage of 
a central heavy ion collision at beam energies Eiab = 5 20 AMeV is studied 
within a microscopic kinetic transport model. The damping rate is shown to 
be a non-monotonic function of Eiab for a given pair of colliding nuclei. This 
fact is interpreted as a manifestation of the zero-to-first sound transition in a 
finite nuclear system. 



I. INTRODUCTION 

The phenomenon of strong dissipation in the collective motions of heated nuclear systems 
is a challenging problem for the nuclear transport theory. Following the analogy with other 
Fermi liquids, like the liquid ^He one expects also in nuclear systems two types of 

collective modes. 

At small temperatures the mean field dominated modes should exist, from the possible 
violatation of the local thermal equilibrium due to the strong Pauli blocking which inhibits 
two-body collisions (the Landau zero sounds). This is actually the main nature of isoscalar 



giant resonances, where indeed the collective energy is essentially given by the amount 



also mean field modes, this time of plasmon type since we can have both a distortion and a 
shift of the neutron (proton) Fermi spheres. 

At high temperatures the usual hydrodynamical collective modes (first sounds) should 
propagate In an infinite Fermi liquid with strong repulsion (Landau parameter Fq ^ 1) 
the two sounds are clearly distinguishable, since the damping rate has a maximum as a 
function of the temperature [03. However, in nuclear matter the transition between the 
two sounds is expected to be smeared-out since the Landau parameter Fq is small in absolute 



The presence of this transition is still an open problem in nuclear dynamics. Only the 
isovector giant dipole resonances can be experimentally studied with sufficient accuracy in 
heated nuclei. The point is that for this kind of two-component mode the transition has 



Aim of this work is to study the transition in the collective nuclear dynamics looking at 
the evolution of the damping mechanism of large amplitude quadrupole oscillations in fusion 
processes, in a microscopic kinetic model. Similar attempts have been recently performed 
in fission dynamics studies We will also follow the temperature dependence of the 

attenuation of the collective mode. 

In the ref.s a reduced friction coefficient 



where is the collective kinetic energy of the dinuclear system (DNS) and {dE / dt) diss 
is the dissipation rate of the total collective energy E = E'f°^ + Epot^ has been extracted 
from the measurements of the prescission neutron multiplicity in fast fission reactions. It 
was shown that in the temperature region T = 2 ^ 3 MeV the reduced friction coefficient 
(5 is very large (/? ~ 10 -i- 100 ■ 10^^ s^^), result not explained with a one-body dissipation 
mechanism only. 



needed to deform the Fermi sphere in momentum space Isovector giant resonances are 



value 0. 



some special features that make it difficult to observe a clear signature [0. 
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In the present work we have performed some kinetic transport studies of the nuclear 
dynamics. Transport equations describe a self-consistent mean field dynamics coupled to 
two-body collisions and so we expect to see in a natural way the transition between the 
two sound propagations. We use the Boltzmann-Nordheim-Vlasov (BNV) ||T0| procedure 
to simulate the phase space dynamics, which has been quite successful in predicting mean 
properties of heavy ion collisions at medium energies. In particular we study central nucleus- 
nucleus collisions leading to fusion at beam energies Eiab < 21 AMeV in order to extract the 
coefficient /5 as a function of the temperature from the damping of the collective quadrupole 
oscillations of the formed di- nuclear system (DNS). Density dependent effective interactions 
of Skyrme type {SKM |[Tl|Jl^) and an averaged free nucleon-nucleon cross section, a = 40m6 
( [|1^) have been used. 

The structure of the work is as follows. In Sect. 2 the procedure of the extraction of the 
reduced damping coefficient (3 from the time dependence of the quadrupole moment given by 
the BNV model is described. In Sect. 3 we present a novel method for determination of the 
temperature based on the energy conservation and on the BNV evolution of the potential 
energy. Sect. 4 contains our results on the /9(T) and their interpretation in terms of two- 
body and one-body dissipation mechanisms. Summary and conclusions are given in Sect. 



II. EXTRACTION OF THE REDUCED FRICTION COEFFICIENT FROM 



According to Eq. (P, the reduced friction coefficient (3 can be, in principle, calculated 
directly, once the phase-space distribution function /(r, p, t) is known from the BNV output, 
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NUCLEAR KINETIC EQUATIONS 



smce: 
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p(r,t) = J dpf{v,p,t) , (4) 
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E,,,{t)=E;:l{t) + - drep{p)p, (5) 



E;il{t)= I dre^.M + E^ui{t) , (6) 



where epip) = ^^/(2m)(37r^p/2)^/^ is the Fermi energy, em.fXp) is the nuclear mean field 
interaction energy density and Ecoui is the Coulomb energy. From Eq.(0) we see that the 
beam energy is not giving contribution to the collective kinetic energy. 

However, in practice, the calculation of the collective kinetic energy E^in is quite am- 
biguous in the test particle technique due to the strong dependence on the width of the 
gaussians representing the test particles . 

We calculate then the coefficient (3 from the time evolution of the quadrupole moment 
of the DNS: 

Qzzit) = Jdr g..(r)p(r, t), g,,(r) = 2^^ - - y\ (7) 

For a damped periodical motion (see Appendix) 

Qzzit) oc exp{—iu!t), u = ujji + iuji, uj < , (8) 

the coefficient /5 is proportional to the imaginary part of the frequency u: 

(3 = -Auji . (9) 

The curves in Fig. 1 show the time evolution of Qzz in the case of central collisions of 
^^Ni + 238^ qx beam energies Eiab = 6.53 20.53 AMeV. The quadrupole moment quickly 
approaches a minimum at t = 100 ^200 fm/c. Afterwards the Qzz starts to grow again, and 
after a time interval At = 50 -h 100 fm/c it saturates. This saturation of the Q^z is caused by 
an almost complete lost of the collective energy of the DNS. One can, therefore, qualitatively 
estimate the [3 coefficient just assuming that \uji\ ur = 2n/tosc, where tosc — 100 fm/c is 
the period of oscillations. That gives f3 = 4:\ujj\ ~ 70- 10^^ s. This value is in agreement with 
the results obtained in Refs. However, we would like to stress that in Refs. ||,^ the 



damping coefficient was extracted from the outgoing part of the trajectory of the DNS from 
compact mononucleus shape to scission during a time interval ~ 30000 fm/c. Our analysis 
is concentrated on the part of the trajectory in vicinity of the mononucleus shape and the 
corresponding time scale ~ 300 fm/c is much shorter. 

In order to get the values of the coefficient /3 we fit a part of the curve Qzzit) to the 
function 

Qif{t)=A + Bsm{ujRt + (l)o)exp{-\iUj\t) . (10) 

The upper time limit of the fitting region is chosen at the second minimum of the quadrupole 
moment. The lower time limit is given by the earliest time when the same value of the Qzz 
as in the second minimum is reached. This definition of the time limits, from one hand, 
corresponds approximately to the selection of a full oscillation in vicinity of the compact 
mononucleus shape. On the other hand, the stage of strong dissipation from the ffist min- 
imum to the ffist maximum is completely in the fitting region. The best fit functions of 
Eq. (|TUp are shown by full dots in Fig.l. The fitting parameters A, B, u^, (pQ, uj and 
the corresponding coefficient (3 = 4:\ujj\ for central collisions of ^^Ni + ^ss^j various beam 
energies are collected in the Table. We see that the damping increases with the increasing 
collision energy until Eiab = 10.53 AMeV. This fact can be understood already from Fig. 1, 
since with increasing beam energy a larger AQzz = — Qmimmum jg (jumped during 

a shorter time interval. But at higher energies Eiab > 10.53 AMeV the damping starts to 
decrease, as we can see from the presence of quadrupole vibrations at later times t > 200 
fm/c (Fig. 1). Eventually at Eiab > 14.53 AMeV a kind of saturation seems to be reached. 

III. DETERMINATION OF THE TEMPERATURE FROM BNV SIMULATIONS 

A direct way to extract the temperature is to fit the local momentum distribution given 
by the BNV model to a T 7^ Fermi distribution. However, in the case of low-energy 
nuclear coUsions studied in present work, this direct method is not appropriate just because 
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its accuracy of ~ 1 MeV is not enough. Therefore we will follow a procedure based on the 
conservation of the energy and on the time evolution of the potential energy given by the 
BNV. The thermal excitation energy of the DNS is (c.f. 

^ therm = -^fcin' ~ ^coul " Erot — ^Epot , (H) 

where -E^j™' = EiabAiA2/ {A1+A2) is the center-of-mass kinetic energy, Ecoui = e^ZiZ2/ {Ri + 
i?2 + 3.5) is the Coulomb energy, = \2A\I^ (fm) z = 1, 2 are the nuclear radii, E,.^, = 
h^L'^ / (26) is the rotational energy with L being the angular momentum (in h units) and 
= [l{AiRl + A2RI) + ^1^2/(^1 + A2){Ri + R2f]mnuc being the momentum of inertia 
{ninuc is the nucleon mass). The last term AEp^t in the r.h.s. of Eq.(|r2) is the difference 
between the values of the potential energy just before and after the overlapping of the density 
profiles. The origin of this term is mostly from the sharp decrease of the nuclear surface 
energy when the two nuclei touch each other. In Fig. 2 we show the time dependence of the 
potential energy per nucleon in the central collision of ^"^Ni + ^38^ 10.53 AMeV. After 
some small increase until t ~ 20 fm/c due to the Coulomb contribution the total potential 
energy quickly drops by about 1 MeV/nucleon reaching the minimum at t ~ 60 fm/c. We 
define AEpot as 

AEp^, = E™" - , (12) 

where E^/^ and -E"™^^ are the first minimum and the first maximum of the potential energy. 

In the Table we report the values of the potential energy "jumps" AEp^, the excitation 
energies and corresponding temperatures T = \J E^f^^^^/ a, a = A/8 MeV~^ for the ^"^Ni + 
^^^U central collisions at various beam energiesQ Our temperatures are higher than those 
obtained in Ref. PJ. In particular, for the Ni+U collision at 6.53 AMeV we have T = 3.1 
MeV and the authors of Ref. |^ report T = 2.4 MeV. This difference is mainly explained 
by the fact, that the temperatures in Ref. [§] are obtained from the outgoing stage of the 



In the case of the central collisions Erot = in Eq. (|Tl]) 
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reaction on much larger time scales as compared to our study. A further contribute to some 
temperature overshooting is coming from pre-equilibrium emissions, neglected in the energy 



balance Eq.(|Il|). This effect, although expected to increase with beam energy, can be still 



considered quite small in this energy range. 

IV. TEMPERATURE DEPENDENCE OF THE DISSIPATION 

Fig. 3 shows the temperature dependence of the reduced friction coefficient (3. The BNV 
results are shown by solid line with dots. The coefficient (3 first increases with temperature 
reaching the maximum at T ~ 4.5 MeV and then at higher temperatures it decreases. 
This signal is indeed quite robust and cannot be related to some overestimation of the 
temperatures at higher beam energies as discussed at the end of the previous section. 

One can interpret the results of the full transport calculations within the approach of 
Ref. 0, based on the analytical solution of the linearized Landau- Vlasov equation 

where Sf and /o are the perturbation and the equilibrium value of the phase space distri- 
bution function, 

SUir, P^t) = J^^J + rmSfir, p', t) (14) 

is the mean field perturbation [p = p/p, p' = p'/p'), g = 4 is the spin-isospin degeneracy of 
a nucleon, A^(0) = gm*pF/(2Tr'^h^) is the level density at zero temperature, Fq and Fi are 
the Landau parameters and m* = m/(l + Fi/3) is the effective mass. 

The collision integral in the r.h.s. of Eq. (|1^) is taken in the relaxation time approxima- 
tion: 

Icoii[5f] ^ -^5/k>2 , (15) 

where 
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SfuM^H E yiUP) I dnp,Yap')5f{v\'=, (16) 

l>2m=-l 

is the part of the perturbation containing the quadrupole and higher multipolarity distortions 
of the Fermi surface. The effective relaxation time r includes two- and one-body dissipation 
contributions: 

The relaxation time T2body was calculated in Ref. for various choices of a nucleon-nucleon 
scattering cross section: 



-1 rr2 



2body 



T'/k (18) 



with K ~ 1900 MeV^fm/c for the isotropic energy independent isospin-averaged nucleon- 
nucleon scattering cross section a^N = 40 mb. For the one-body relaxation time we have 
used the wall-and- window formula (c.f. Ref. p): 

'^Ibody — ' ^^^^ 

where R = 1.2A^^^ (A = Ai + A2) is the radius of a mononucleus composed of the two 
colliding nuclei, 
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~^ 6 \eF 



(20) 



is an average velocity of nucleons, and is a numerical factor which depends on the multi- 
polarity and on the isospin of a collective mode. We have chosen a value ^ = 1.85, which 
corresponds to the isoscalar quadrupole mode in the scaling wall model (see Ref. and 
refs. therein). 

The solution of Eq. (p!3| ) inside a nucleus with uniform nonperturbed density can be 
found clS db superposition of plane waves. In this CclSG, clS it was shown in Ref. ||15|| , the 
intrinsic width of a giant multipole resonance F = 2h/Trei, where Trei is the relaxation time 
of the distribution function 6f {6f oc exp(— icut), u = Ufi + iuj, uj = — l/r^ez), can be 
expressed as 

8 



with q = [|(H-Fo)(H-Fi/3)]~^. Eq. describes both well known zero sound (T oc r^^, 
ljrt ^ 1) and first sound (F oc r, cu/jr ^ 1) regimes. In our calculations we have put the 
Landau parameters Fq = 0.2, Fi = 0. The frequency has been chosen as Hur = 64.7^4"^/^ 
MeV {A = Ai + A2) corresponding to the giant quadrupole resonance. 
The solid line in Fig. 3 shows the friction coefficient 

PzFST = 2T/h (22) 



with F given by Eq. (|2l|). We see that Eqs. (ID agree qualitatively with the BNV 



model, but the analytical calculation gives a more smeared-out transition between the two 
sounds. 

In the limit of the zero sound the coefficient (3 is 

Pzs = 4/r = (3ffy + Pl'fy , (23) 

where jS'^g'^^ = 4:/T2body, Pz's'^^ — ^/'^Uody This simple formula (dot-dashed line in Fig. 3) 
is quite close to the BNV results at T < 4.5 MeV. However, either the two-body (dashed 
line) or one-body (dotted line) contributions taken separately are strongly underpredicting 
the BNV friction coefficient. 

In order to better understand the relative importance of the two-body and one-body 
mechanisms we have also performed the BNV calculations at Eiab = 6.53, 8.53, 10.53, 12.53 
and 14.53 AMeV switching off the collision term, corresponding to a pure Vlasov evolution. 
Fig. 4 shows the comparison of the Qzz time evolution with and without collision term 
for the Ni+U collision at 14.53 AMeV. In the case without collisions the damping is much 
reduced: we observe several large- amplitude oscillations of the quadrupole moment. We 
have also checked that in the Vlasov evolution the oscillations at later times, t > 200 fm/c, 
are present for all the other studied beam energies. In Fig. 3, the coefficient (3 in the 
case of coUisionless dynamics is shown by full squares. The damping slowly increases with 
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temperature and saturates at T = 5 MeV. As expected, the system is always in the region 
of the zero sound. The absolute value of the reduced friction coefficient in the collisionless 
case is 3 -^ 5 times less than for a full BNV calculation and close to the wall-and-window 
value 

V. SUMMARY AND CONCLUSIONS 

The reduced friction coefficient was extracted for the initial stage (~ 300 fm/c) of central 
^^Ni + 238u collisions at beam energies Eiab = 6^21 AMeV in the framework of the BNV 
transport model. The quadrupole moment Qzz has been chosen as a relevant collective 
variable. The quadrupole time evolution shows overdamped oscillations with a damping 
rate proportional to the friction coefficient. Two-body collisions play a major role in the 
damping of Q^z- 

As a function of the beam energy, the damping rate has a maximum at Eiab — 10 AMeV. 
The corresponding temperature of the DNS, neglecting particle emissions, is 4.6 MeV. We 
interpret, therefore, this temperature as a transition temperature from the zero-to-first sound 
propagation. 

This result agrees with earlier calculations of V.M. Kolomietz et al. (Ref. |^), where 
a value of the transition temperature 4^5 MeV was deduced on the basis of the analyt- 
ical solution of the linearized BUU equation in the case of isoscalar giant resonances of 
multipolarities / = and 2 in a hot nucleus. 

Our calculations are overpredicting the transition temperature 2 2.5 MeV that follows 
from the analysis of the prescission neutron multiplicities by J. Wilczyhski et al. (Ref. 
0). We have to stress that the considered collective modes are different, in our analysis 
quadrupole oscillations in the entrance channel dynamics while in the ref. ^ the fission 
mode in the exit channel. In particular a relevant variance is on the time scales of the 
two modes, with a much larger proper time for the fission dynamics. Following the simple 
condition ujrt ~ 1 for the transition, we roughly get a Ttr oc ^/uJR and then we can expect 
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quite smaller transition temperatures for the fission modes. Moreover, as already stressed 
before, likely our temperature assignements are a little overestimated, particularly for the 
higher beam energies, due to the lack of pre-equilibrium emissions in the energy balance. 

Looking at the BNV simulations, the transition from zero-to-first sound appears as a 
presence of quadrupole vibrations at relatively late times, t > 200 fm/c (see Fig. 1), for 
beam energies above 10 AMeV. These vibrations would accelerate the fission of the produced 
DNS that should be observed experimentally, as an increase of fast-fission cross sections. 
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APPENDIX 



In this Appendix we derive the relation between the imaginary part uj of the complex 
frequency to = + iuji and the reduced damping coefficient (3 of Eq. (|I]). 
From the Eq. and the continuity equation p = — V(vp) we get: 



Qzz = J py^Vq,,d'r . (24) 
In the case of the small- amplitude damped periodical motion of the kind 

v(r, t) = (5vo(r) exp(— ztut) (25) 
the Eq. (p^) can be linearized with respect to small values of 6vq: 



Qz 



Po J 6vo{r)Vqzzd^r 



exp{—iujt) = Qzz(t = 0) exp{—iujt) , (26) 



where po is the nonperturbed density (p = Po + ^p)- The general solution of Eq. (^) is 

Qzz{t) = const + exp{-iujt) . (27) 

UJ 
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Therefore in a linear approximation, over Svq, the quadrupole moment reveals oscillations 
with the same frequency uj = ujr + iuji as the velocity field v(r, t). 
The time derivative of the collective kinetic energy (0) is: 



Ekin ^ mpo j Re(v)Re(v)rf^r 



mpo 



5vl{v)d?r 



cos{ujRt){uJi cos{ujFit) — ujr sin(co'/jt)) exp(2c<j/t) . (28) 

Averaging Eq. ( ^8]) over the period of oscillations we come to the relation: 

Ekin ^ ujipom J (Re(v))2c?V = 2uJi'Ei~ , (29) 

where we have dropped the therm oc cos{ujjit) sin(co'^t) exp(2co'/t) which changes the sign 
during the period of oscillation. According to the virial theorem for harmonic oscillators 

Ekin = Epot = -E , (30) 

where E is the total collective energy. Therefore 

E ~ 4u;jE^ . (31) 
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Table Fit parameters from Eq. (|TOD; potential energy "jumps" AEpot] excitation energies 
and temperatures of the DNS Eqs. (|1T]) , (|12|) ; for central collisions ^^Ni + ^38^ 
various beam energies Eiab- 



Elab 

(AMeV) 


A 
(Afm2) 


B 

(Afm2) 


c/fm 


00 


c/fm 


-AEpot 
AMeV 


E* 
AMeV 


T 
MeV 


6.53 


33.81 


24.34 


0.024 


-6.537 


0.0133 


0.88 


1.2 


3.1 


8.53 


25.60 


39.94 


0.030 


-6.169 


0.0175 


1.28 


1.9 


3.9 


10.53 


21.46 


38.45 


0.036 


-6.197 


0.0180 


1.68 


2.7 


4.6 


12.53 


18.34 


31.13 


0.038 


-6.091 


0.0163 


2.10 


3.4 


5.2 


14.53 


16.42 


27.03 


0.039 


-5.906 


0.0147 


2.28 


3.9 


5.6 


16.53 


14.33 


23.50 


0.040 


-5.843 


0.0137 


2.46 


4.4 


6.0 


18.53 


12.89 


26.73 


0.039 


-5.589 


0.0141 


2.68 


5.0 


6.3 


20.53 


11.14 


28.05 


0.038 


-5.486 


0.0147 


2.74 


5.4 


6.6 
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FIGURE CAPTIONS 



Fig. 1 Solid lines - time dependence of the quadrupole moment Qzz for ^^Ni + ^38^ central 
collisions at beam energies (from top to bottom) Eiah = 6.53, 8.53, 10.53, 12.53, 14.53, 
16.53 and 20.53 AMeV. The value of the quadrupole moment at t=0 is always 85 
fm^/nucleon. Full circles show the best fits obtained using Eq. ( [T(]| ) for each collision 
energy. See text after Eq. (PUj ) for a definition of the fit regions. 

Fig. 2 Potential interaction energy per nucleon Eq. as a function of time for the 
64Ni(10.53 AMeV) + ^^^U central colhsion. 

Fig. 3 Temperature dependence of the reduced friction coefficient (3. BNV calculations 
with (without) collision term are shown by full circles (squares) connected with thin 
solid line to guide eye. Errorbars are due to the ambiguity caused by a finite number 
of test particles. The analytical result of Ref. for the Giant Quadrupole Resonance 
in a hot nucleus (see Eqs. (pT]), (p^)) is shown by thick solid line. The two- and 
one-body contributions and their sum (see Eq. (^3])) in the zero sound limit are shown 
by dashed, dotted and dash-dotted lines respectively. 

Fig. 4 Comparison of the time evolution of the quadrupole moment for ^^Ni(14.53 AMeV) 
_j_ 238^ central collision calculated with (solid line) and without (dashed line) collision 
term. 
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